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Abstract Nonlocal currents, in devices where two 
normal-metal terminals arc contacted to a superconduc- 
tor, are determined using the circuit theory of meso- 
scopic superconductivity. We calculate the conductance 
associated with crossed Andreev reflection and electron 
transfer between the two normal-metal terminals, in ad- 
dition to the conductance from direct Andreev reflection 
and quasiparticlc tunneling. Dcphasing and proximity 
effect are taken into account. 



1 Introduction 

Transport between a normal-metal and a superconduc- 
tor at subgap energy is possible through Andreev reflec- 
tions, where an incident electron from the normal-metal 
is retro-reflected as a hole and a Cooper pair is trans- 
ferred into the superconductor pQ. However, since An- 
dreev reflection is a nonlocal process on the scale of the 
coherence length, the retro-reflected hole can end up in 
another normal-metal contacted to the superconductor. 
This process is known as crossed Andreev reflection (CA) 
[2[3] and contributes to a nonlocal conductance. We de- 
fine the nonlocal conductance in a three-terminal device 
(see Fig. Q]) as the current response in one normal-metal 
terminal (Ni) to a voltage bias between another normal- 
metal terminal (N2) and a superconducting terminal (S). 
The nonlocal conductance will also depend on electron 
transfer (ET) between the two normal-metal terminals, 
which gives a contribution with opposite sign as CA to 
the nonlocal conductance!]] The nonlocal conductance of 
such systems has recently been studied extensively both 
experimentally [MHZ] and theoretically [5MITO1J1T1JT21 
l4l ll3[fT4l[l~5] since it demonstrates inherently mesoscopic 



physics, and since crossed Andreev reflection is a way to 
produce spatially separated entangled electron pairs. 




1 In our previous paper Ref. [4] this transport process was 
referred to as electron cotunneling, but as this phrase applies 
to the tunneling limit we will here use the more general term 
electron transfer. 



Fig. 1 Our circuit theory model: A cavity (c) is connected to 
one superconducting (S) and two normal-metal terminals (Ni 
and N2). The three connectors are described by their sets of 
transmission probabilities. A coupling to ground represents 
the "leakage current" (see text). 



In our recent paper Ref. [4], we used the circuit the- 
ory of mesoscopic superconductivity [16] to calculate 
the conductances in a three terminal device where two 
normal-metal terminals and one superconducting termi- 
nal are connected to a region where chaotic scattering 
takes place. We assumed that that the energy of the in- 
jected particles was much smaller than the gap of the 
superconducting terminal. In this case, the only trans- 
port process involving only one of the normal-metal ter- 
minals and the superconductor is direct Andreev reflec- 
tion (DA), where the hole is backreflected into the same 
normal-metal as the incident electron. We now extend 
this approach to take into account situations where the 
bias voltage is comparable to the gap of the superconduc- 
tor so that incident electrons from a normal-metal may 
be transferred into the superconductor as quasiparticles 
(QP). Taking these processes into account, the current 
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at energy E out of Ni can be written 

h{E ) = [i _ fl[E ) - f 2 (-E)] 



Get{E) 



[h(E) - h(E)] 



e 

Gqp(E) 



[f s (E) - h(E)] 



(1) 



Here, the functions f n (±E) denote the Fermi-Dirac dis- 
tribution functions in terminals n = 1,2, S at energy 
±E and we have defined energy dependent conductances 
G{E) for the transport processes discussed above. Total 
charge current is given by /charge, l = J dEIi(E). The 
nonlocal differential conductance is obtained from Eq. 

o, 

<9^chare;e,l 



= - j AE [G ET (E) - G CA (E)} 



df{E-eV 2 ) 



dE 



(2) 



This shows how the nonlocal conductance is determined 
by competing contributions from crossed Andreev reflec- 
tions and electron transfer. 



2 Circuit theory 

In the circuit theory of mesoscopic superconductivity 
[16] . a system can be modeled as a network of termi- 
nals, connectors and nodes in a similar manner as elec- 
tric circuits on the macroscale are modeled using classi- 
cal circuit theory. The theory is formulated in terms of 
matrix Green's functions and matrix currents in Nambu- 
Keldysh spac^l and e.g. describes quantum effects and 
flow of charge, energy and particle-hole correlations. 

To describe the three-terminal devices at hand we 
introduce equilibrium Green's functions depending on 
temperature and local chemical potential for the two 
normal-metal terminals and the superconducting termi- 
nal. The terminals are coupled to a region where scatter- 
ing takes place, which we will refer to as a cavity. Our 
circuit theory model is shown in Fig. [1] If we assume 
that the distance between the contacts is very small on 
the scale of the coherence length, the properties of the 
cavity are spatially homogeneous so that it can be de- 
scribed by one nonequilibrium Green's function G c . If 
spatial variation in the scattering region is important, it 
can be modeled as a network of cavities with different 
Green's functions. We assume that G c is isotropic due 
to chaotic or diffusive scattering. 

Connectors between terminals and the cavity are de- 
scribed by their sets of transmission probabilities {T^}, 



In our notation a% is a Pauli matrix in Nambu space and 
fi is a Pauli matrix in Keldysh space. Matrices with structure 
in both spaces are denoted with accent ". 



with flow of matrix currents /„ depending on the Green's 
functions of adjacent elements, 



k 



(n) 



[G„ , G c 



4 + Tf ) ({G„,G c }-2) 



(3) 



The matrix current describes not only the flow of charge, 
spin, and energy currents, but e.g. also the flow of cor- 
relations [To] . 

Finally, the theory is completed by a generalized "Kirch- 
hoff 's rule" : The sum of matrix currents flowing into a 
cavity should vanish. This determines the nonequilib- 
rium Green's function of the cavity in our present circuit. 
The spectral charge current through connector n can 
be determined 7x, n = Tr j<7 3 i r l K ' j /8e once the Green's 
functions G„, G c have been determined. Similarly, the 
spectral energy current becomes /l,™ = Tr|ii K '| /8e. 

The K superscript denotes the Keldysh matrix block of 
the current. 

Using the generic circuit theory model above, we can 
describe crossed Andreev reflection in a wide range of 
different experimental systems. For example, consider 
a system where two adjacent normal-metal electrodes 
are deposited on and connected by metallic contacts 
to a mesoscopic superconductor. The terminal Green's 
functions correspond to properties of the normal-metals 
and the superconducting sample away from the contact 
region. The cavity corresponds to the small region of 
the superconductor between the normal-metal contacts. 
Metallic contacts are modeled by putting Tjf 1 ' = 1 for 
the conducting modes of contacts n = 1,2. The contact 
between the cavity and the superconducting terminal 
can be modeled as a diffusive connector introducing a 
bimodal distribution of transmission probabilities [17j . 
Other systems for experimental study of crossed An- 
dreev reflection may be fabricated from superconductors 
coupled to semiconductors where the geometry is defined 
by deposition of gates. For example, the semiconductor 
can consist of a ballistic cavity with two point contacts 
to normal reservoirs. In this case, the cavity Green's 
function describes the nonequilibrium state of the bal- 
listic cavity. The point contacts are modeled by putting 

(n) 

Tf = 1 for the open channels and zero otherwise for 
n = 1,2. The transparency and number of conducting 
modes of the contact to the superconductor can also be 
determined experimentally. 

Superconducting pairing and dephasing in the cavity 
are described in circuit theory by introducing a "leakage 
current" /leakage = ^ ie 2 i / oV c [G c , H c ] in the matrix cur- 
rent conservation on the cavity [16j . Here i>q is the den- 
sity of states, V c the volume, H c = Ea^ + i<7iRc {A c } + 
io^Im {A c } the Hamiltonian, and A c the gap in the cav- 
ity. The energy-dependent term in H c describes dephas- 
ing between electrons and holes, and sets an energy scale 
for the proximity effect which we will refer to as the ef- 
fective Thouless energy of the cavity. When we disregard 
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Josephson effects, the phase of A c can be chosen arbi- 
trarily e.g. purely imaginary and inspection of the re- 
tarded part of the matrix current conservation reveals 
that pairing inside the cavity appears with the same 
matrix structure as the coupling to the superconduct- 
ing terminal. Therefore, the effect of pairing inside the 
cavity can be described by a renormalization of the cou- 
pling between the cavity and the superconducting ter- 
minal. Thus the difference between a normal and a su- 
perconducting cavity is equivalent to rescaling this con- 
ductance, i.e. only quantitative modifications. 



processes: 



3 Results 



We will now discuss the structure of the Green's func- 
tions and the solution of the matrix equations. The nor- 
mal terminals have Green's functions Gin) — + 
(o"3/iL,i(2) + 1^t,i(2))(t = i +1^2), where we have introduced 
the charge- and energy-distribution functions /it {E) , /il {E) 
that can be written in terms of the particle distribu- 
tion function f(E) as /it = 1 — f{E) — f(—E) and 
h L = -f(E) + f(-E), sec Rcf. [18]. The retarded (ad- 
vanced) part of the Green's function of the supercon- 

([E±iS}as+Aia2) / O , where 



ducting reservoir is G^ A ^ 



Q = y/{E ± id) 2 - A 2 , A is the gap. The Keldysh part is 
obtained from Gf- = Gg hg — hgGg where h$ = ct 3 /il,s + 
I/iT,s- We parametrize the Green's function of the cavity 
as Sf = (T 3 cosh(6») + i<7 2 sinh(0) and G£ = -<7 3 (G*)t<7 3 , 
the Keldysh part is given by Gf = G* h c - h, ^ A 

h c = thh.c + (T3/lT,c- 



G^, where 



With the Green's functions specified as above, we im- 
pose matrix current conservation in the cavity, ^2 n I n + 
^leakage = to obtain equations that determine 0, h^ c , 
and /il,c- The equations for the distribution functions 
are conservation of charge (energy) at each energy, 



(L),n (,"T(L),n 



*T(L) 



(4) 



for n = 1,2, S, where we have defined effective, energy 
dependent conductances for charge (energy) G'Y^- )n [8(E)] 
between reservoir n and the cavity. Using the solution for 
G c that we have obtained, we calculate the current out of 
terminal Ni and compare the result to Eq. (Q]). This al- 
lows us to determine the conductances for the transport 



Gqp(E) 
Gba(E) 



G ET (E) 

CA 



Gl,iG 



L.S 



Gl,1 + Gl,2 + Gl,S 
1 / Gt,1 (Gt,2 + Gt,s) 
4 



Gt,1 + Gt,2 + Gt,S 
Gl,1 (Gl,2 + Gl,s) 
Gl,1 + Gl,2 + Gl,S 
1 ( _ _ Gl,iGl,2 
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± 



Gl,1 + Gl,2 + Gl,S 

Gt,iGt,2 



G 



T,l 



r T,2 



G 



T,S 



(5a) 



(5b) 



(5c) 



The conductance for quasiparticlc transport into the 
superconducting terminal, Gqp, is proportional to Gl,s 
which is the energy conductance for transport from the 
cavity into the superconductor. This quantity vanishes 
at subgap energies where Andreev reflection of particles 
from opposite sides of the Fermi surface is the only pos- 
sible transport process. The symmetry between Get and 
Gca in Eq. (|5c|) was discussed in our previous paper Ref. 
|4), and shows that the differential nonlocal conductance, 
given by Get — Gca see Eq. ([2]), is always positive. In 
that paper, we also discussed the limit that there is no re- 
sistance between the superconducting reservoir and the 
cavity. We see from Eq. (|5c|) that in this case Get(ca) 
vanishes because of the large terms Gt,s> Gl,s- 

The equations that determine 9 must generally be 
solved numerically. We have performed such calculations, 
and show in Fig. [5] the result for a system where connec- 
tors to Ni and N 2 are point contacts and the connector 
to S has transparency ; = 0.5 for the conducting 
modes. Defining g n = e- 2 ^2 k T^ / {-kK) we choose pa- 
rameters gj/gs = 0.1 for j = 1,2. There are now two 
energy scales in the problem, A and -Erh, and we have 
chosen A/ETh = 6 in Fig. [21 The nonlocal conductance 
is largest for energy smaller than the effective Thoulcss 
energy, defined as i?Th = %s/(2e 2 i'oV r c ), and decreases 
in two steps at £rh and A with increasing energy. The 
nonlocal conductance above the gap corresponds to the 
normal state result dlxjdVi = g\gil{g\ + .92 + ,9s) for 
quasiparticlc transport in a three terminal network. In 
this energy range we also have a contribution from quasi- 
particlc transfer into S. 

When the transmission probabilities of the interface 
to S are in the tunneling limit, i.e., all C 1 the CA 
conductance will be suppressed in comparison to the case 
shown in Fig. [2] The ET conductance, on the other hand, 
is enhanced by the reduced transmission of the contact to 
S. In Fig. [3] we show the conductances for a device where 
Ni and N 2 are connected by point contacts, and S by 
a tunneling barrier. The conductance of the tunneling 
barrier is the same as the conductance of the contact 
to S in Fig. [21 For energies between £>rh and A, we 
see that Gca is strongly suppressed. The total nonlocal 
conductance, Get — Gca, has a minimum in the subgap 
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Fig. 2 Conductances when Ni and N2 are connected by 
point contacts, and S by an interface of intermediate trans- 
parency Tj_ — 0.5. The nonlocal conductance is largest for 
energy smaller than the Thouless energy, and has features 
at energy corresponding to the gap of the superconducting 
reservoir at A = 6 2?Th- 



regime at energy corresponding to £rh, and is largest 
for an energy close to A. This is qualitatively different 
from the device where S is connected by an interface of 
intermediate transparency (Fig. [2]), where the maximum 
nonlocal conductance in the subgap regime is found at 
very small energy, and then decreases with increasing 
energy due to increasing CA conductance. 
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Fig. 3 Conductances when Ni and N2 are connected by 
point contacts, and S by a tunnel barrier. In contrast to the 
the device where S was connected by an interface of interme- 
diate transparency (Fig. [2}, the CA conductance below the 
gap is strongly suppressed in this case. 



4 Conclusion 

In conclusion, we have studied nonlocal transport in a 
three-terminal device with two normal-metal terminals 



and one superconducting terminal. To this end we have 
applied the circuit theory of mesoscopic transport. The 
connectors between the circuit elements are represented 
by general expressions, relevant for a wide range of con- 
tacts. Dephasing is taken into account, and gives rise 
to an effective Thouless energy. We calculate the con- 
ductance for crossed Andrcev reflection, electron trans- 
fer between the normal-metal terminals, and direct An- 
drcev reflection and quasiparticle transport between one 
normal-metal terminal and the superconducting termi- 
nal. The nonlocal conductance is generally dominated 
by electron transfer in our model, similar predictions 
were made in Refs. (12114"] . We showed in Ref. [I] that 
for this model, in the limit that there is no resistance 
between the device and the superconducting terminal, 
our results agree with Ref. [8] and the total nonlocal con- 
ductance vanishes. We numerically calculate the conduc- 
tances for experimentally relevant combinations of con- 
tacts to demonstrate the appearance of two energy scales 
in the conductances: The effective Thouless energy and 
the gap of the superconducting terminal. The conduc- 
tance for crossed Andreev reflection depends strongly 
on the transparency of the interface to superconducting 
terminal as demonstrated in our numerical calculations. 
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